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We search for the rare leptonic decay B + — » t + u t in a sample of 232 x 10 6 BB pairs collected 
with the BABAR detector at the SLAC PEP-II B-Factory. Signal events are selected by examining 
the properties of the B meson recoiling against the semileptonic decay B~ — > D*°£~Vi. We find no 
evidence for a signal and set an upper limit on the branching fraction of B(B + — > t + v t ) < 2.8 x 10 -4 
at the 90% confidence level. We combine this result with a previous, statistically independent BABAR 
search for B + — » t + v t to give an upper limit of B(B+ -> t+v t ) < 2.6 x 10" 4 at the 90% confidence 
level. 

PACS numbers: 13.20.He, 14.40.Nd, 14.60.Fg 



In the Standard Model (SM) the purely leptonic decay 
B + — > t + v t 1] proceeds via the annihilation of the b and 
u quarks into a virtual W boson. Its amplitude is propor- 
tional to the product of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix [2| element \V u b\ and the B meson decay 
constant fs ■ The SM branching fraction is given by Q : 

B(B+ t+ Vt ) = ^^ml (l - f%\V ah \*T B , 

(1) 

where Gp is the Fermi coupling constant, m T and mj 
are the r + lepton and B + meson masses, and tb is the 
B + lifetime. The branching fractions for B + — > e + v e 
and B + — * p^v^ are helicity-suppressed by mf/m B , 
where mg is the mass of e + or Using the value of 
\V ub \ = (3.67 ±0.47) x 10" 3 and the lattice QCD cal- 
culation of f B = (0.196±0.032) GeV 0, we determine an 
expected value of B(B+ -> t+v t ) = (9.3 ± 3.9) x 10~ 5 . 
Currently, our best knowledge of fg comes from theoret- 
ical calculations, with a current theoretical uncertainty 
of roughly 16% J3]. Observation of B + — ► t + v t could 
provide the first direct measurement of fs ■ The ratio of 
B(B + — > t + v t ) and Am^, the difference in heavy and 
light neutral Bd masses |6| , can be used to determine the 
ratio of CKM matrix elements jKi&|/|Vtd| with roughly 
4% theoretical uncertainties [4J, |i| , dominated by the un- 
certainties on square root of the bag parameter \fB~B [E| . 

No evidence of the B + — ► t + v t decay has been re- 
ported to date. The most stringent published experi- 
mental limit is B{B+ -> t + v t ) < 4.2 x 1(T 4 at the 90% 
confidence level (C.L.) [jj. Physics beyond the SM, such 
as supersymmetry or two-Higgs-doublet models, could 
enhance B(B + — > t + v t ) up to the current experimental 
limits 0. 

The data used in this analysis were collected with the 
BABAR detector || at the PEP-II asymmetric-energy 
e + e~ storage ring. The results are based on a data 
sample of (231.8 ± 2.6) x 10 6 BB events, in an inte- 
grated luminosity of 210.6 fb^ 1 collected at the T{AS) 
resonance. An additional sample of 21.6 fb _1 was col- 
lected at a center-of-mass (CM) energy approximately 
40 MeV below the T(45) resonance. We used the latter 
sample to study continuum events, e + e~ — > qq (q = u, d, 



s, c) and e + e~ — ► t + t~. Charged-particle tracking and 
dE/dx measurements for particle identification (PID) 
are provided by a five-layer double-sided silicon vertex 
tracker and a 40-layer drift chamber operated in the 
1.5 T magnetic field of a superconducting solenoid. A 
detector of internally reflected Cherenkov light (DIRC) 
is used to identify charged kaons and pions. The energies 
of neutral particles are measured by an electromagnetic 
calorimeter (EMC) consisting of 6580 CsI(Tl) crystals. 
The magnetic flux return of the solenoid is instrumented 
with resistive plate chambers in order to provide muon 
identification. A full detector Monte Carlo (MC) sim- 
ulation based on EvtGen [irl | and GEANT4 [nj is used 
to evaluate signal efficiencies and to identify and study 
background sources. Beam-related background and de- 
tector noise samples are obtained from random triggers 
at regular intervals. These samples are overlaid on the 
simulated events with appropriate luminosity weighting 
to model these time- varying background conditions. 

Due to the presence of at least two neutrinos in the 
final state, the B + — ► t + v t decay lacks the kinematic 
constraints that are usually exploited in B decay searches 
in order to reject both continuum and BB backgrounds. 
The strategy adopted to search for this decay is to recon- 
struct the B - meson from an T(AS) — ► B + B~ event in 
a semileptonic final state, denoted by B~j". All remain- 
ing charged and neutral particles in that event, referred 
to as the "signal-side" particles throughout this paper, 
are then examined under the assumption that they are 
attributable to the decay of the accompanying B + ("sig- 
nal B"). 

The B^ is reconstructed in the decay modes B^ — > 
D*°£~V e {£ = e or //). The D*° is reconstructed in 
the modes B 7r° and D ^. The D° is reconstructed in 
four decay modes: K~ir + , X~7r + 7r~7r + , K~7r + ir°, and 
K®ir + TT~. All kinematic variables are calculated in the 
CM-frame of the Y(4S) unless otherwise noted. 

Photon candidates are obtained from EMC clusters 
with laboratory-frame energy E 1 greater than 30 MeV 
and no associated charged track. Photon pairs with in- 
variant mass between 115 and 150 MeV/c 2 are taken as 
7T° candidates. 
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The D° candidates are reconstructed by selecting com- 
binations of identified pions and kaons with invariant 
mass within 40 MeV/c 2 of the nominal D° mass ex- 
cept for the K~Tt + Tr° mode, where this window is 70 
McV/c 2 . Each D° candidate is combined with a soft ir° 
or 7 candidate to form a D*°. The ir° and 7 candidates 
are required to have momentum less than 450 MeV/c. 
Further, the 7 candidate must have E y > 100 MeV. The 
invariant mass difference AM between the D*° and D° 
is required to be within the range 135-150 MeV/c for 
the D°tt° mode, and 130-155 MeV/c 2 for the D°j mode. 

The £Tj" — > D*°£~17£ candidates are identified by com- 
bining a D*° candidate of momentum p D ,« > 0.5 GeV/c 
with a lepton candidate of momentum pi > 1.0 GeV/c. 
The lepton candidate must be identified as either an elec- 
tron or a muon. The invariant mass mu-oi of the D*°£ 
candidate is required to be greater than 3.0 GeV/c 2 . Un- 
der the assumption that a massless neutrino is the only 
missing particle, the cosine of the angle between the di- 
rections of the -BT and the lepton-£>*° combination is 

2 Sboam • Eo*°i — m B — m 2 D , p 

co S 6 BtD , H = — . , , (2) 

2\p D ,o e \ ■ V^ hcam -m B 

where -Ebeam is the expected B~ meson energy. The en- 
ergy and momentum of the D*°£ candidate are E D *og 
and Po*<>ei respectively. Correctly reconstructed candi- 
dates populate the range [—1, 1], whereas combinatorial 
backgrounds can take unphysical values well outside this 
range. We retain B~\ candidates in the wider interval 
|cos 0B.D*°e\ < 1-1, allowing for the effects of detector en- 
ergy and momentum resolutions. If more than one D*°£ 
candidate is reconstructed in an event, the best candi- 
date is selected using a likelihood based on the simulated 
D° mass and AM distributions. We further require that 
the sum of the charges of all the particles in the event 
( "net charge" ) must be equal to zero. 

The B~, reconstruction efficiency for events contain- 
ing a B + — > t + v t decay is determined from signal 
simulation after verifying that the simulated BB, uu, 
dd, ss, cc, and t + t~ events are consistent with data. 
This procedure compensates for differences in the 
reconstruction efficiency in the low-multiplicity environ- 
ment of B + — > t + v t events compared with the generic 
B + B~ environment. The simulated efficiency is fur- 
ther cross-checked by comparing the yield of events in 
which a B + — + D* £ + vi decay has been reconstructed 
in addition to a B^ ("double semileptonic decays"). In 
the signal simulation the B^ reconstruction efficiency 
is £ s i = (1.75 ± 0.07(stat.) ± 0.05(syst.)) x 10~ 3 . The 
D* £~17i, D*°, and D° branching fractions are factored 
in £ 8l . 

Events that contain a -ETj" are examined for evidence of 
a B + — > t + v t decay. Charged tracks and EMC clusters 
not already utilized for the B~ x reconstruction are as- 
sumed to originate from the signal candidate B + decay. 



We identify the r lepton in six mutually exclusive chan- 
nels: e + v e T> T , p + v^U Tl 7r + j7 T , 7r + 7r°I/ r , 7r + 7r~7r + i7 r , and 
"misidentified lepton" . The misidentified-lepton channel 
selects signal events from the e + VeP T or /i+j^^v signal 
decays in which the momentum of the e + or p + from the 
signal t + is too low to pass the lepton identification cri- 
teria. The identified t + modes all together correspond 
to approximately 81% of all t + decays 0. 

Signal candidates are searched in events that are re- 
quired to possess exactly one signal-side charged track, 
except for 7r + 7r~7r + 77 r candidate events, which must have 
three signal-side charged tracks. The signal track from 
the e Jr v e y T {p + v^V T ) channel is required to be identified 
as an electron (a muon), and not to satisfy either muon 
(electron) or kaon PID criteria. In the 7r + 77 r , 7r + 7r°F T , 
7r + 7r~7r + 77 T , and misidentified-lepton channels the signal 
track(s) must not satisfy electron, muon, or kaon PID. In 
addition, each signal track from the 7r + 7r~7r + I7 r channel 
has to be identified as a pion. For the 7r + 7r°77 r chan- 
nel the signal track is combined with a signal-side ir° 
candidate, reconstructed from a signal-side photon pair 
(E~ ( > 50 MeV for each photon) with invariant mass 
between 100 and 160 MeV/c 2 . If several signal-side ir° 
candidates are reconstructed in an event, the candidate 
with 77 invariant mass closest to the nominal ir° mass Q 
is chosen. We require that the events in the 7r + 77 r and 
misidentified-lepton channels contain no signal-side ir° 
candidates. Events in the Tt + v T and misidentified-lepton 
channels are distinguished by requiring the momentum 
of the signal track to be greater than 1.2 GeV/c in the 
former, and less than 1.2 GeV/c in the latter. 

Further requirements are made on the (total) momen- 
tum of the signal track(s) for some channels: p e + < 
1.4 GeV/c for e + v e V T , and p n + n -^+ > 1.0 GeV/c for 
7r + 7r~7r + 77 T . We apply constraints on the missing mass 
M m i ss of the event, which is determined by subtract- 
ing the total four-momentum of reconstructed tracks 
and neutrals from that for the Y(AS) system. This 
quantity tends to be larger for events with more neu- 
trinos. Signal events must satisfy M m i SS > 4 GeV/c 2 
for e + is e 17 T and /i + ^I7 T , M m i ss > 3 GeV/c 2 for n + V T , 
TT + Tr°17 T and misidentified-lepton, and M m j ss > 2 GeV/c 2 

for 7T + 7T~7r + F T . 

Additional kinematic constraints are applied on the 
ir + ir°v T (■k + ~k~'k + T> t ) channel, which proceeds mainly via 
intermediate p + (af and p°) resonance(s). In the 7r + 7r°F T 
channel the invariant mass of the ir + ir° must be between 
0.55 and 1.0 GeV/c 2 . For the 7r + 7r~7r + 77 r channel the 
invariant mass of the three-pion system is required to be 
within the range 1.0-1.6 GeV/c 2 . The ir + Tr~ combination 
of the three-pion system, with invariant mass closest to 
the nominal p° mass 4], is required to have momentum 
greater then 0.5 GeV/c and invariant mass between 0.55 
and 1.0 GeV/c 2 . We further require that the cosine of 
the angle between the directions of the r + and the 7r + 7r 
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(ir + TT 7T + ), 

COS T , had = 



_ 2E T ■ E had - ml- ml ad 
2|Pt| ■ |Phad| 



(3) 



is within [—1.1, 1.1]. Here -Ehad, Phad and Whad are the 
energy, momentum and invariant mass, respectively, of 
the 7r + 7T° (7r + 7r~7r + ). The energy E T and momentum p r 
of the r + from B + — > r + ^ r decay are calculated under 
the assumption that the B + is at rest in the CM-frame. 

Continuum background events contribute to the 7r + F r , 
misidentified-lepton, ir + ir°17 T , and n + n~n + v T channels. 
To suppress this background we combine five variables 
in a linear Fisher discriminant [l2|: Pd'°, Pi, co s®B,D'°ti 
the cosine of the angle between the thrust axis of the 
decay products of B~ { and the thrust axis of the rest 
of the event, and the ratio of the second and zeroth Fox- 
Wolfram moments using all the particles in the event . 
The requirement placed on the output of the Fisher dis- 
criminant selects about 93% of signal events and rejects 
about 37% of continuum background events. After this 
requirement the continuum background in each channel 
is less than 40% of the total background. 

The sum of the laboratory-frame energies of the neu- 
tral EMC clusters with E 1 > 30 MeV, which are not 
associated with either the £Tj" or the ir° candidate from 
7r + 7r°F T channel, is denoted by E extIa (Fig.HJ. For signal 
events the neutral clusters contributing to B e xtra come 
only from hadronic shower fragments, bremsstrahlung, 
and beam-related background. This variable peaks near 
zero for signal while for background, which contains ad- 
ditional sources of neutral clusters, it takes on larger 
values. Signal events are required to have E extra less 
than 250 MeV for e+v e V T , 150 MeV for n+VpVr, 300 
MeV for n + V T , 170 MeV for misidentified lepton, 250 
MeV for n + Tr°V T , and 200 MeV for 7r+7r"7r+I7 r , which 
are selected based on MC study to provide the tightest 
branching fraction upper limit. The -Eextra selection re- 
gion defines the "signal region" for each channel. The 
350 < E cxtla < 1000 MeV region is defined as the "side 
band" for all the channels. 

The efficiencies £j for each r selection channel i are de- 
termined using simulated events. Cross-feeds among the 
r-decay channels are taken into account. The system- 
atic uncertainties in the selection efficiency arise from 
tracking efficiency (1.4% per track), particle identifica- 
tion (0.2%-2.0%), Bextra simulation (3.0%-8.0%), tt° re- 
construction (3.3%), and data and MC differences in the 
output of the Fisher discriminant (1.0%). Systematic un- 
certainties due to the -Eextra simulation are determined 
by evaluating the effect of varying the MC -E C xtra distri- 
bution within a range representing the observed level of 
agreement with data in samples containing B~y and up 
to seven additional tracks. For further cross-check the 
-Eextra distributions of the data and MC events for the 
double semileptonic decays are compared. The signal se- 
lection efficiencies for the six selection channels are listed 
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FIG. 1: The distribution of __/ ex tra after applying all other se- 
lection criteria, plotted for (a) e + v e u T , (b) /i + i/ M 77 T , (c) 7r + 77 T , 
(d) misidentified lepton, (e) 7r + 7T 77 T , and (f) 7r + 7r tv + v t 
channels. The data and background MC samples are repre- 
sented by the points with error bars and solid histograms, 
respectively. The dotted lines indicate the B + — > t + v t 
signal distribution from MC. The signal MC events for the 
e + v e T T , n+VfipT, 7T + 77 T , and misidentified-lepton (tt + 7v°17 t 
and 7r + 7r~7r + !7 r ) channels are normalized assuming a branch- 
ing fraction of 10 -3 (10 -2 ) for B + — > t + v t decay. 



in Table [I] The total B + — > t + v t selection efficiency is 
roughly 31%. 

The remaining background consists primarily of B + B~ 
events with correctly reconstructed Z?T . For these events 
the signal side contains fC°'(s), neutrino(s), or particles 
that pass outside the detector acceptance. For each chan- 
nel we estimate the background 6j in the signal region us- 
ing events in the data side band and the simulated E ext r& 
distribution: 



Ardata (pj-MC /atMC \ 



(4) 



Here N^^b ^ s the number of data events in the side band, 
and A^gf^ and NgfdeB are the numbers of MC background 
events in the signal region and side band, respectively. 
Background estimation is cross-checked using data and 
MC events that satisfy the full signal selection, with the 
exception of having two signal-side tracks, or non-zero 
net charge, or the AM of the D*° outside the selection 
region. The uncertainties in the background estimations 
are predominantly statistical; smaller systematic uncer- 
tainties arise from the simulation of the -E C xtra shape in 
the background MC. 

We determine the B + — > t + v t branching fraction from 
the number of signal candidates Si expected for each r 
selection channel, where Sj = N B ± e s \Ei B(B + — > t + v t ). 
N B ± = (231.8 ± 2.6) x 10 6 is the estimated number of 
S ± mesons in the data sample. The results for each 
channel are combined using the estimator Q = C(s + 
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TABLE I: Efficiency (si) with statistical and systematic er- 
rors, expected background (bi), and observed data candidates 
(rii) for each reconstructed r selection channels. The cross- 
feeds among the r decay modes are taken into account. The Si 
values include the branching fractions of the r decay modes. 



selection 



£,(%) 



7.5 ±0.4 ±0.2 13.4 ±2.4 17 

2.9 ±0.2 ±0.1 6.2 ±1.7 5 

8.0 ±0.4 ± 0.3 27.7 ±5.0 26 

2.5 ±0.2 ±0.1 28.6 ±4.3 31 

1.4 ±0.2 ±0.1 21.6 ±3.0 26 

misidentified lepton 9.0 ± 0.4 ± 0.4 33.4 ± 5.1 45 



7T + 7r ■K + V T 



b)/C(b) HIHH, where 



db' 



(5) 

is the likelihood function for signal-plus-background hy- 
potheses, rii is the observed number of data events in each 
t selection channel, and di is the uncertainty in the back- 
ground estimate bi (Table QJ. The likelihood function for 
background-only hypotheses C(b) can be obtained from 
Eq. El by setting Si to zero. 

The measured branching fraction, which is the 
value that maximizes the likelihood ratio estimator, is 
(1.3i"*'i) x 10~ 4 - This value is compatible with a zero 
branching fraction. The and 6, values (Table HJ do 
not indicate any significant excess of observed events. 
Therefore, we set an upper limit on the branching frac- 
tion [3 of B{B+ -> t+v t ) < 2.8 x 1CT 4 (90% C.L.). The 
expected branching fraction upper limit for background 
only hypothesis is B{B+ -► t+v t ) < 1.8 x 1(T 4 (90% 
C.L.). 

The BABAR Collaboration has previously performed a 
search for the B + — > t + v t decay based on a sample of 
88.9 x 10 6 BB pairs, where the B~ meson accompanying 
the signal B + is reconstructed in a variety of hadronic 
or semileptonic modes Q. The hadronic B~ selection is 
mutually exclusive with the current B~^ selection. There- 
fore the two samples are statistically independent and 
may be combined. The hadronic reconstruction analy- 
sis obtained a limit B(B + — ► t + v t ) < 4.2 x 10~ 4 at 
the 90% C.L. To combine the results from the previous 
hadronic and current semileptonic samples, we create a 
combined estimator from the product of the semileptonic 
(Qsi) and hadronic (Qhad) likelihood ratio estimators, 
Q = Qs\ x Qhad- The measured branching fraction from 
the combined sample is (1.3l 'g) x 10~ 4 . This value is 
compatible with a zero branching fraction, and we set a 



combined upper limit, 

B(B+ -> t+v t ) < 2.6 x 10~ 4 (90% C.L.). (6) 
These results represent the most stringent limits on 
B + — ► t + v t reported to date. 
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